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ABSTRACT 



Simultaneous measurements of sea surface elevation and onshore and 
alongshore water particle velocities were measured at three locations 
within the surf zone using two capacitance type pentrating wave staffs 
and three two-component electromagnetic flow meters. The probability 
density functions, pdf, for the sea surface elevation were always highly 
positively skewed, whereas the pdf's for the velocities were both nega- 
tively and positively skewed. Mean values of the onshore and alongshore 
components of flow reflected the influence of a rip current frequently 
observed just south of the instrument locations. Strong harmonics in 
the spectra of sea surface fluctuations and particle velocities infer 
nonlinear conditions. Coherence values between waves and onshore flow 
were high, ranging above 0.9. The coherence between waves and onshore 
flow was used to separate the turbulence and wave-induced velocity 
components. Over the range of collapsing to spilling breakers a rea- 
sonable value for the ratio of turbulent to wave-induced velocity was 
determined to be approximately 0.75. Saturation regions were found in 
the wave and velocity energy-density spectra at higher frequencies as 
evidenced by -5 and -3 slopes, respectively. 
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I. INTRODUCTION 



A. HISTORICAL PERSPECTIVE 

The surf zone is an area bounded on the seaward side by the point 
where waves first begin to break, and on the landward side by the point 
of maximum run-up on the beach slope. Although beach profiles and waves 
within the surf zone have been studied by many, the kinematics of wave 
forms and water particle velocities in the zone have remained something 
of a mystery. 

The problems encountered have been practical as well as theoretical. 
Wave theories developed for deep water waves do not correctly characterize 
the motion which occurs as the waves break. Direct measurements of break- 
ing waves appear to be the most viable means of approaching the problem. 
However, this has been hampered by the surf zone's hostile environment, 
which is difficult to reproduce in the laboratory, and by inadequate 
instrument design. Improvements in the latter area has resulted in sturdy, 
sensitive, but expensive measuring devices which have rapid response times. 

The earliest study of significant importance was conducted by 
Iversen (1953) in which he used photographic techniques to obtain a 
Lagrangian description of water particle motion under breaking waves. 

The laboratory channel limited the wave type to plunging and surging 
breakers . 

Inman (1956) was one of the first to make sophisticated field 
measurements, measuring the drag force on a cantilevered sphere to infer 
water particle motion. Miller and Zeigler (1964) measured in situ 
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particle motions using acoustic and electromagnetic current meters, 
then compared their findings with higher order wave theory and found 
some qualitative agreement. Walker (1969) made similar measurements 
using propeller type flow meters. Meadows (1976) used ducted impeller 
flow meters at equally spaced vertical positions in order to measure 
longshore components of flow. Huntley (1976) utilized a single two- 
component electromagnetic flow meter in an effort to obtain a value 
of the friction coefficient. Bowen and Huntley (1974) made measure- 
ments of the nearshore velocity fields using as many as three two- 
component electromagnetic flow meters. Fiihrbbter and Biisching (1974) 
measured simultaneous orbital velocities and water levels using a two- 
component current meter and two pressure type wave meters. Thornton 
(1969), Steer (1972), Thornton and Richardson (1973), Bub (1974) and 
Galvin (1975) used pressure meters, capacitace wave gauges and electro- 
magnetic current meters tocompute particle velocities and surface profiles 
within the surf zone. The work presented herein is an extension of these 
latter studies. 

B. OBJECTIVE 

The objective of this research is to study the kinematics of water 
particles within breaking waves and within the surf zone. In particular, 
the experiment was performed on a flat beach in order to measure the 
characteristics of spilling breakers. Simultaneous measurements were 
made of the instaneous sea surface elevation and two orthogonal water 
particle velocities, onshore and alongshore, at three fixed locations 
in the surf zone. Estimates of the probability density functions and 
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power spectra of the wave heights and particle velocities are made. The 
suitability of using linear wave theory as a spectral transfer function 
in determining velocity spectral components from the power spectrum of 
the waves is measured. Computed velocity spectra are compared with 
actually measured velocity spectra. Coherence between waves and onshore 
flow is used to separate the turbulence and wave-induced velocity com- 
ponents in order to determine the ratio of turbulent to wave-induced 
velocity. Spectral estimates are analyzed to ascertain the slope of 
the saturation region at higher frequencies. 



11 



II. MEASUREMENTS 



A. EXPERIMENT 

The experimental site is just north of La Jolla, California. Data 
presented in this study was taken as part of a larger scale experiment 
conducted from 15 February to 30 March 1977 in the Southern California 
bight area. The over-all purpose of the experiment was to test and 
evaluate the instrumentation package which was selected for use aboard 
SEASAT. The apparatus, which will measure waves, winds, and sea surface 
temperature, was being flown on board NASA airplanes. A major ground 
truth program was concurrently being conducted at Torrey Pines Beach. 

Measurements of waves, currents, and sediment transport at Torrey 
Pines Beach were made using the following equipment. A five element 
linear array of pressure sensors was positioned offshore along the ten 
meter contour to measure wave energy and direction. Located in the surf 
zone inshore of the array were three support towers in a line perpendi- 
cular to the beach. A surface piercing capacitance wave gauge and a two- 
component electromagnetic flow meter were mounted on each tower. A 
nephelometer was installed on the middle support tower. Tc the north 
of the three towers was an array of four flow meters which was used to 
delineate the areal and temporal variation in alongshore currents. 

Figure 1 shows instrument locations. Resistance-wire run-up meters were 
installed in the swash zone at four locations. All instrumentation 
was battery operated. 

The offshore pressure sensor outputs were telemetered directly to 
the shore processes lab at Scripps Institute of Oceanography , which is 



12 



15ft- 



0 



O 

oo 

on 



oo 

OT 



0 



o 

uo 



cc 

o 



00 

CD 



O 

<c 



oo 



CD 



<C 

or 



o 

cn 



oo 

or 



o 

z 

c 

oo 

UJ 

CD 

=) 

c 

CD 






o 

LO 



00 



OO o 

h- 

c nr 

DD O 
cd nr 



4r 



4- 



(SU313W) 3U0HSdd0 39NVH 



c o 
>- 
q s: 
nr <c 

>— « CJ 



0 ■ 



13 



Figure 1. Instrument Location for Overall Experiment. 



located one mile to the south. All data from the other instruments 
were cabled to one of two transmitting terminals on the beach where it was 
then telemetered to the shore processes lab and recorded. All the equip- 
ment except the nephelometer were in position prior to high tide on 9 
March. The nepholometer was installed on 14 and 15 March. Data was 
accumulated during high tide on 9, 10, 11, 16, 17, 18, 19, 21 and 23 
March. 

Several other smaller scale experiments were conducted at the same 
time. At an elevation of approximately 300 feet on a cliff overlooking 
Torrey Pines Beach a radar installation was established. The radar 
images of waves were then photographed and analyzed. 

Sediment transport experiments were conducted on 11, 21, and 23 
March. Bed load transport was measured using the methods described by 
Inman and Komar (1970) to trace the movement of fluorescent dyed sand. 
Suspended sediments were measured in situ by swimmers using a mechanical 
water sampling devide. These samples were taken along with nepholometer 
readings in order to compare methods. During these sediment transport 
investigations, Lagrangian floats made from wine bottles weighted by 
sand were used to determine average long-shore current speeds. 

B. INSTRUMENTATION. 

The manner in which waves break depends very much on the character- 
istics of the beach and near shore bottom slope. (Table I) Spilling breakers 
occurred most frequently at the Torrey Pines Beach site. Figure 2 shows 
a typical beach profile and the location of the instruments which provided 
the data for this paper. 
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TABLE I. BEACH AND WAVE CHARACTERISTICS 



LOCATION 

DATE 

BEACH SLOPE 
SAND TYPE 
SAND MEDIAN GRAIN 
BREAKER TYPE 
WAVE PERIOD 
WAVE HEIGHT 



TORREY PINES BEACH 
9-23 MARCH 1977 
FIGURE 2 and APPENDIX A 
QUARTZ 

SIZE 1.2 MILLIMETERS 

SPILLING 

11.9 to 15.9 SECONDS 
UP TO 1.5 METERS 
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FIGURE 2 . Typical Beach Profile and Instrument Location at Torrey Pines Beach 

on 19 March 1977 



1 . Wave Gauge 



The wave gauges were of the capacitance type fashioned from 
3/8 inch outside diameter stainless steel rod. The rod was tightly 
covered with 1/16 inch wall thickness polypropyl ene tubing. This li- 
near, highly sensitive instrument has proven to be sufficiently sturdy 
to withstand the rigors of the unfriendly surf zone. The gauges operate 
on the principle that a change in the plate dimension of the capacitor 
changes its capacitance and consequently the circuitry output voltage. 

In these gauges, the insulated steel rod and sea water act as the plates 
and the insulation functions as the dielectric. As the water level fluc- 
tuated, the capacitance changed. Fluctuations were sensed by a transis- 
torized circuit powered from the beach. The circuit was designed by 
McGoldrick (1969). The electronics packaged was housed in a watertight 
brass case which was mounted on the tower. This allowed the connecting 
leads to be relatively short, hence minimizing wire-to-wire capacitance. 
All gauges were statically calibrated in the laboratory prior to the ex- 
periment. Accuracy was estimated to be -.005m. The calibration plots 
are shown in Appendix B. 
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2. Flow Meter 



The flow meters were Marsh-McBirney model 511 electromagnetic 
water current meters. Operation of the meter is based on Faraday's prin- 
ciple of electromagnetic induction. The meter measures water particle 
velocity in two orthogonal directions through a range of zero to three 
m/sec with a maximum output error of two percent of full scale reading. 
The sensor has a variable time constant, a setting of 0.2 sec. was used. 

The flow meters were dynamically calibrated using the method of 
Thornton and Krapohl (1974). This procedure involves utilizing an os- 
cillating platform driven by a variable speed motor geared to an eccen- 
tric throw arm. This technique was employed in order to determine the 
meter's characteristics. Measurement accuracy was determined to be 
*.005 m/sec. The only difficulty encountered in their utilization oc- 
curred when Flow Meter #1 was buried on 19 March due to the heavy sedi- 
ment transport. 

The measuring instruments were attached to a 3.6m high tower con- 
structed of steel pipe with an outside diameter of 6. 3cm. A 0.5m diameter 
base plate was placed about 0.6m from the bottom of the tower. This 
configuration allowed the tower to be sunk approximately 0.7m into the 
sand. The towers were also supported by steel guy wires fastened to 
blade anchors driven into the sand. Thus, tower movement and vibration 
were negligible. 

The towers were placed on a line perpendicular to the shore at 
specified distances, and were erected during low tide when the beach 
was more accessible. As indicated earlier, measurements were taken at 
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high tide. Instruments were arranged so as to be in the same vertical 
plane. A carpenter's level was used to establish axis alignment with an 
estimated error of -2 degrees. A typical tower and instrument arrangement 
is shown in Figure 3. 



19 



CAPACITANCE 
WAVE GAUGE 




20 



FIGURE 3. Schematic of Typical Instrument Tower with Sensors 



III. ANALYSIS OF DATA 



The data were digitized at a rate of four samples per second, 
corresponding to a sample interval of 0.25 seconds. This resulted in a 
Nyquist frequence of 2.0 Hz. The sampling rate was considered suffi- 
ciently high enough to avoid aliasing of energy into the portion of the 
spectra which was of interest. 

Continuous time series records of all measurements were taken for 
roughly 2-1/2 hours each day. Record lengths of 24 minutes from each 
data set were subjected to analysis. Criteria which determined record 
lengths included economy of computer usage and resolution over the fre- 
quency range of interest. 

A maximum lag time was chosen as 5 % of the record length. This re- 
sulted in a spectral band width resolution of 0.007 Hz and each spectral 
estimate having 40 degrees of freedom. The ninety percent confidence 
limits for 40 degrees of freedom using a chi-square distribution are found 
to be between 0.72 and 1.51 of the measured power spectral estimates. 

A mean value was computed for all data sets and the data was linearly 
detrended to remove tidal effects. Variance, standard deviation and 
average period were calculated. The average period was determined by 
calculating the time between zero upcrossings. The calculated average 
periods are lower than the visually observed periods due to a number of 
perturbations, such as noise and capillary waves, increasing the zero 
upcrossing occurrences. A probability density function for each data set 
was calculated and plotted. Comparisons were made with Gaussian and Gram- 
Charlier distributions using the chi-square goodness-of-fi t test. Variance, 
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standard deviation, skewness and kurtosis of the distributions were 
computed and are summarized in Table II. 

For each data set an auto-covariance function was determined and then 
smoothed with a Parzen window. Applying a Fourier transform to the smoothed 
auto-covariance function resulted in the power spectrum. An examination of 
the power spectrum shows the regions of greatest potential or kinetic 
energies and their respective frequencies. 

Cross-spectra were computed by calculating a cross-covariance between 
data sets, smoothing it with a Parzen window and applying a Fourier trans- 
form. Coherence and phase were then determined from the cross-spectrum. 

A coherence and phase versus frequency plot will indicate the regions 
and degree of linear relationship and phase between two data sets. 
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TABLE II. COMPUTED STATISTICS 
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TABLE II. COMPUTED STATISTICS (Continued) 

STANDARD CHI-SQUARE PARAMETER 

DATE INSTRUMENT MEAN VARIANCE DEVIATION SKEWNESS KURTOSIS GAUSSIAN GRAM-CHARLIER 

T9 - March WG#1 1.25410 0.03791 0.19472 1.15469 5.92272 2524.97 1431.11 
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Mean and variance for the wave heights is given in meters squared; mean and variance for 
the velocity components is given in meters squared per second squared. 

Standard deviation for the wave heights is given in meters; standard deviation for the 
velocity components is given in meters per second. 



IV. RESULTS 



A. QUALITATIVE DESCRIPTION 

Observation of various breaker types which occur on different beaches 
has resulted in a number of similarities being noted. A typical analog 
record of waves and onshore velocities obtained at Torrey Pines Beach is 
shown in Figure 4. In general, a sawtooth shaped profile is created when 
there is a quick drawdown of water just before the breaker arrives, fol- 
lowed by a steep, vertical leading edge, and a sloping profile toward 
the trailing edge. On the trailing edge, secondary waves are often ob- 
served. Secondary waves are harmonics of the primary wave frequency and 
denote strongly nonlinear waves. At Torrey Pines Beach the breaking 
waves are generally of the spilling variety which spill rapidly at the 
crest and move down the wave. 

B. MEAN VALUES 

A mean value was computed for all data sets and linearly detrended 
to remove tidal effects. The sea surface elevation mean values ranged 
from 0.67 to 1.33. Wave heights increased till the 18th and then de- 
creased during the remaining period data was taken, reflecting well the 
actual environment. 

Analysis of the mean values for both the onshore and alongshore 
components of flow was complicated by an apparent rip current. A rip 
current was frequently observed located just south of the array of wave 
gauges and flow meters. The on-offshore flow for the five days under 
consideration at Flow Meter #1 was always directed offshore at about 0.1 
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ALONGSHORE VELOCITY 

Figure 4. Typical Analog Record of Waves and 
Velocities from Torrey Pines Beach. 
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m/sec. At Flow Meter #2 the horizontal water particle velocity was di- 
rected onshore at roughly 0.25 m/sec. The flow at Flow Meter #3 was 
away from the coast on 16, 17, and 18 March at slow speeds and towards 
the shore on 21 and 23 March moving very slowly. Concurrently , the long- 
shore flow at Flow Meters #1, #2 and #3 was generally moving downcoast 
with the peak velocity varying from flow meter to flow meter from day to 
day. Similar results were reported by Huntley and Bowen (1974) for a 
nearshore circulation cell which they attributed to a rip current. 

C. PROBABILITY DENSITY FUNCTIONS 

1 . - Gaussian and Gram-Charl ier Frequency Distributions 

The computed probability density function was compared to the 
Gaussian and Gram-Charl ier distributions and tested for comparison using 
the chi-square goodness-of-fit test. The closer the chi-square fit 
parameter is to zero, the better the fit (Table II). 

The Gaussian probability density function is given by 

P G a (x ) = 1 e -(x-u)W (,) 

a /2 tT 

where a is the standard deviation and y is the mean. The Gaussian 
distribution is completely described by the mean and variance, and has 
zero skewness and a kurtosis equal to three. This approximates the 
sea surface in deep water conditions. However, nonlinearities found in 
the surf zone introduce skewness and kurtosis values that deviate from 
Gaussian and result in a distribution more closely approximated by the 
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Gram-Charl ier distribution. This is to be expected in that the Gram- 
Charlier pdf is calculated using the assigned parameters from the data, 
and includes moments higher than the second. The frequency distribution 
of Wave Gauge #1 on 23 March is shown in Figure 5. Appendix D contains 
distributions calculated for other data sets. 

Distributions for all flow velocities are either unimodal or have 
less pronounced secondary peaks compared to the distributions of wave 
heights. The distribution of the onshore flow of Flow Meter #1 on 23 
March is shown in Figure 6. As noted in Table II, this results in nearly 
equal fit parameters with the Gram-Charl ier distribution being generally 
smaller and thus a better fit. Fit parameters for velocity components 
are much closer to zero than are the values calculated for wave components 
since the flow meters do not experience the surface irregularities to 
the same degree as do the wave gauges. 

2. Skewness and Kurtosis 

The sea surface distribution was found to be positively skewed. 

This means that positive values are large, but less frequent than negative 
values which are more frequent, but smaller, indicating a greater amount 
of time below the mean water level. In other words, the crests are steeper 
and more peaked, whereas the troughs are elongated and flatter which is 
how waves in shallow water theoretically appear. 

It should follow that the wave-induced particle velocities would 
have distribution similar to the waves. However, of the 34 cases listed 
in Table II, 17 had positive skewness and 17 had negative skewness with no 
readily apparent pattern. Skewness of the onshore flow at Flow Meter #1 



28 



GRAM-CHARLlLR DH-TR1BUT10 



r^ 

cn 




cj 

s*_ 

ro 



co 

cvi 



c 

o 



=H= 

CU 

cn 

3 

ro 

<L3 

CD 

> 

ro 



C 

o 



ro 

> 

CU 



CU 

u 

ro 

<*- 

s- 

3 

UO 

ro 

CU 

GO 

CU 



c 

o 



3 

.o 

s- 

+-> 

cn 

O 

>> 

C_) 

c 

CU 

3 

O' 

cu 



<u 

i- 

3 

cn 




30 



-002 -nm ooo ooi 002 003 

STANDARD DEVIATIONS 

Figure 6. Frequency Distribution of the onshore flow at Flow Meter #1 on 23 March 1977. 



was always positive. Alongshore flow at the same location had negative 
skewness in three of four cases. Onshore water particle velocity at 
Flow Meter #2 had negative skewness for all six days. The skewness of 
alongshore flow at the same point was positive in five of six cases. 

The onshore component of flow at Flow Meter #3 had a positive skewness 
for 16, 17, and 18 March and a negative skewness for the following three 
days on which data was taken. Skewness of the alongshore flow at Flow 
Meter #3 was positive for four of the six days. The reason for this 
disparity is not known. Similar anomalous results have been found by 
Thornton and Galvin (1975). 

Kurtosis indicates the peakedness of each parameter. Just prior 
to breaking the waves achieve the greatest degrees of peakedness. Visual 
observations indicated that most of the waves were breaking at or near 
Wave Gauge #1 on 17, 18, and 19 March; this correlates with their high 
kurtosis values on these days. Values of kurtosis for all flow veloci- 
ties again had no apparent pattern with values ranging from 2.6 to 8.6. 
The reason for this is not known. 

D. SPECTRAL ANALYSIS 

1 . Sea Surface Elevation and Water Particle Velocity 

The power, coherence, and phase spectra of Wave Gauge #1 and the 
onshore component of the flow of Flow Meter #1 on 23 March are shown in 
Figure 7. The spectrum of the wave surface elevation, which is char- 
acteristic of all the records, show a narrow banded peak at a frequency 
of 0.07 Hz, corresponding to a period of 14.3 sec. Also evident are 
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Figure 7. Power, Coherence and Phase Spectra for the onshore flow of 
Flow Meter #1 and Wave Gauge #1 on 23 March 
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peaks at 0.14 and 0.21 Hz which appear to be harmonics of the primary 
peak at 0.07 Hz. These harmonics appear to be physical as was observed 
from the strip chart earlier. The harmonics probably have some energy 
contribution due to the Fourier computational technique. Appendix E 
contains spectra calculated for all other wave and onshore flow data sets. 

The coherence values were high, ranging to greater than 0.9 in 
the maximum energy portion of the wave. Coherence begins to fall off 
at about 0.5 Hz except on the 17th and 18th of March, the two days on 
which short-crested waves were noted. On these days the coherence begins 
to decrease at about 0.3 Hz. The decrease in coherence is probably due 
to the relative increase in turbulence, the general decrease in energy 
level, and wave energy spreading. The decrease in coherence due to wave 
energy spreading is discussed in Section E. 

The phase angle, which according to linear theory should be zero, 
had an average value of less than 20 degrees over the highly coherent 
band of prominent wave energy. The phase angle tends to increase slightly 
over the coherent band and becomes random for the noncoherent region of 
the spectrum. 

In order to compare waves and flow velocity, the wave profile 
spectrum was converted to a theoretical velocity spectrum for comparison 
with the measured flow velocity spectrum using the linear transfer function 
such that 

y f ) = | H ( f ) | 2 S n (f) (2) 

The transfer function is given from linear wave theory 
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H(f) = 2 irf coshk(h+z) 
sinh kh 



( 3 ) 



where 

k is the wave number, 

h is the mean water depth, and 

z is the depth of flow meter below the mean depth (m) 

Quantitatively, the results showed that linear theory overestimated 
wave-induced velocity spectral components by about 50%; in previous 
studies by Thornton et_ al_ (1976) where the waves were long crested and 
arrived at near normal incidence, the opposite was true, in that linear 
theory underestimated the magnitude of the wave spectrum. The difference 
in results is probably due to wave directionality, since when the wave 
approaches at any angle other than normal the flow meter measures only a 
component of the actual flow. 

2. Onshore and Alongshore Flow 

The cross-spectra between the onshore and alongshore components 
of flow were calculated with the intent of finding the relation between these 
two parameters. Power, coherence, and phase spectra of the onshore and 
alongshore flows at Flow Meter #3 on 23 March are shown in Figure 8. 

The spectra, which are representative of all the records, shows very 
little coherence between velocity components. The phase angle is ran- 
dom over the noncoherent band. 

The reason for the lack of coherence between the two orthogonal 
flow components is believed to be due to wave directional ity . As was 
seen earlier the onshore component of flow appears highly wave-induced; 
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Fiqure 8 Power, Coherence and Phase Spectra for onshore and 
alongshore flow at Flow Meter #3 on 23 March. 
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concomitantly, the alongshore flow component is at near normal incidence 
to the wave direction and thus has only a small wave-induced contribution. 
The variance of the alongshore component of f-low is significantly less than 
the on-offshore component except at zero frequency, which reflects the 
mean longshore current variation. 



E. TURBULENT VELOCITY VERSUS WAVE-INDUCED VELOCITY 



The total velocity can be separated into components of a mean, plus 
wave-induced, plus turbulent velocity 

u = u + U + u' . (4) 

The wave-induced and turbulent velocity spectral components are assumed 
to be statistically independent. For unidirectional waves aligned with 
the measured velocity component, the co-spectra between waves and velocity 
is then given by 

V f > = W f >- < 5 > 

Assuming statistical independence, the horizontal velocity spectrum 
is obtained from 



s u( f ) = S u ,(f) + S u (f). (6) 

A further assumption is made that the waves and wave-induced velocities 
are given by a constant parameter linear process where the coherence 
is equal to unity 



*U f) 



5 U7?(f) 



= 1 . 



(7) 



V f )V f) 

Substituting (5), (6), and (7) into the definition of coherence between 
the total horizontal velocity and waves results in 
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-1 




i + V(f) 

S U (f) 




s 



U 



( 8 ) 



Increasing lack of coherence is due to an increasing high ratio of turbu- 
lence. The coherence values indicate the percent of total velocity 
which is associated with the wave velocity. 

Wave-induced velocity can be determined using (8) 



The wave-induced velocity spectral component calculated in this way 
should be less than the actual value, the reason for this being that 
the linear coherence between waves and velocities will always be under- 
predicted due to nonlinearities which are always large in a breaking 
wave and also as a result of directional spreading of the incident 
wave energy (Battjes (1974) ). 

The turbulent velocity spectrum is derived by subtracting the cal- 
culated velocity spectrum from the measured velocity spectrum. 



The ratio of the turbulent velocity intensity to the wave-induced 
velocities intensity is shown in Figure 9. The data used for collapsing 
and plunging breakers, denoted by squares, was measured during earlier 
studies by Thornton et al_ (1976). The collapsing and plunging breakers 
were measured on beaches where the waves were long crested and arrived 
at near normal incidence. The stars represent data collected during this 
experiment. A value of the velocity intensity is given by taking the 



V f > r, (f)s u (f). 



(9) 



s u ,(f) - S u (f) - s 0 (f). 



( 10 ) 
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Figure 9. Ratio of turbulent to wave-induced velocity intensities 



square root of the variance which is computed by integrating the 
velocity spectrum, S u (f), across all frequency bands 

CT u'C^i< f > df - 01) 

The measurements show a small spread of values for the ratio of turbulent 
to wave-induced velocity intensities. The average ratio calculated is 
approximately 0.82. This would indicate that, to a first approximation, 
the wave-induced kinetic energy is approximately 1.2 times the turbulent 
kinetic energy. Thus the velocities in the surf zone on the average are 
highly wave-induced. 

If waves are not unidirectional, and generally they are not, it is 
necessary to account for directional spreading. The wave angle striking 
the wave staff and the type of wave, long or short crested, must be con- 
sidered. In general, short crested waves predominated at Torrey Pines 
Beach. On the 17th and 18th of March they were particularly noticeable. 

The coherence between surface elevation and velocity components of short 
crested waves is less than that of long crested waves due to directional 
spreading. Using the wave train model developed by Yefimov and Khristoforov 
(1971) the ratio of coherence between waves and onshore flow for the short 
crested waves to that for long crested waves is given by 

V^f.Os) (1 + l/3ctn 2 0 o + 4S u ,(f)/3S u (f) ) _1 

4 f .01 ) ( 1 + S .(fJ/SylfJCl/cos! ] r 1 

' 0 

The model describes long crested waves using a train of two-dimensional, 

planar random waves traveling in direction 9 0 ■ This corresponds to a 

bivariate spectrum given by 
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(13) 



i|> n (f,01) = S n (f)6(0-0o) , 



where 6(0-0 o ) is the Dirac delta function. 

Short crested waves appear as a train of three-dimensional waves 

whose average direction of spectral components is assumed to be of the form 

( 2/tt ) S (f) cos 2 (0-0 o ) , 

0 for| 0 - 0 O |> it /2 , (14) 

Neglecting turbulence in order to examine the effect of directionality, 
the ratio of the coherence between waves and wave-induced velocities for 
short crested waves to that of long crested waves becomes 



(f»0 s ) 



yfal f,0 s ) (1 + l/3ctn 2 0o) _1 

y lb ue ' > 1 



(15) 



As an example, when do equals an angle of 30°, the ratio equals 0.9. 



Conversely, the influence of turbulence on the ratio of coherence 
between waves and velocities for short crested waves to that of long 
crested waves can be examined by disregarding the directionality and 
letting 0 o =O. The ratio is given by 



ry.fls) = 1 + 4S u ,(f)/3S u (f) )' 

^( f ,0i) ( i + s u .(f )/s u (f ) )■ 



-1 



( 16 ) 



A summary of the ratios considering the turbulence and/or direction- 
ality factors at various angles is tabulated in Table III. The turbulence 
factor, S y ,Cf )/Sy ( f ) , is assumed for comparison purposes to be equal to 
0.75, the average value for the measured data in Figure 9 as denoted by 
the squares. 
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Table III 
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2 2 
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M f * ,/ y 4 f i> 


y „„( f .fl s )/y JfA) 

u>7 s U !j i 




(Equation 15) 


(Equation 16) 


(Equation 12) 


0 




0.875 


0.875 


5 


0.997 




0.877 


10 


0.990 




0.882 


15 


0.977 




0.891 


20 


0.957 




0.905 


25 


0.932 




0.923 


30 


0.9 




0.947 



The coherence between waves and velocities was used to calculate the 

wave-induced velocity (equation 11). This method under-predicts the actual 

wave-induced velocity because of nonlinearities and wave energy spreading. 

Since the turbulent velocity is computed by subtracting the calculated 

wave-induced velocity spectrum from the measured velocity spectrum 

(equation 12) the smaller wave-induced velocity value makes the turbulent 

velocity value larger. The average ratio of the turbulent velocity 

intensity to the wave-induced velocity intensity was calculated to be 

0.825. In order to correct for the directional spreading it is necessary 

2 2 

to multiply by the ratio 7 u ^(f, 0 S )/V (f,#-|). As an example, the 

ratio of the coherences for an angle of 15° is equal to 0.89 1 . Multiplying 
this value by the average ratio results in a turbulence to wave-induced 
ratio of 0.75, which is equal to the earlier measured mean value. Hence 
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over the range of collapsing to spilling breakers a reasonable value for 
cr 1 /cr ( j is approximately 0.75. 

F. SATURATION REGION IN THE SPECTRUM OF BREAKING WAVES 

Waves on Torrey Pines Beach characteristically broke as spilling 
breakers with the crest sliding down the face of the wave. When the 
fluid particles at the free surface move forward at a greater speed 
than the wave speed, c, breaking occurs. Strong harmonics observed 
in the wave and velocity spectra (Figure 7) are indicative of energy 
being transferred from low to higher frequencies. The energy is even- 
tually dissipated by viscosity at the highest frequencies. When the 
transfer of energy is not rapid enough to balance the increase in energy 
density of the waves during shoaling, breaking occurs. The waves are 
said to be "saturated" with energy. Hence, a region of saturation 
would be expected through which energy is transferred from the low to 
higher frequencies. During shoaling and breaking of waves on a beach, 
the waves become saturated at low frequencies first imposing a bound 
on the peak energy density. 

Spectra of waves at breaking, or just inshore of breaking, are 
shown in Figures 10 and 11 plotted on a log-log scale. The spectral 
estimates have been normalized by dividing by the total variance in an 
effort to bring the high end of the spectra to a single line. The 
spectral estimates have been block-averaged over frequency bands such 
that the log-energy values are linearly distributed on the log-frequency 
axis. The confidence intervals decrease with increasing frequency 
because of the increasing number of frequency bands averaged over. The 
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confidence interval is noted on the figures. The slopes of the log-log 
spectra at high frequencies for Wave Gauges #1 and Wave Gauges #3 can 
be approximated by a -5 slope. This is indicative of a saturation or 
equilibrium region for wind generated waves in deep water. This was a 
surprise, as a -3 slope was expected for the waves measured in shallow 
water, Thornton (1976). A possible explanation is that most of the waves 
measured were waves inside the initial break point and would be classified 
as spilling breakers or reformed waves. It may be that these types of 
waves do not reach saturation condition as found for the more intense 
plunging and collapsing breakers. 

The normalized velocity spectra are given in Figures 12 and 13. The 
slopes of the velocity spectra at higher frequencies at Flow Meter #1 and 
Flow Meter #3 are both closely approximated by a -3 slope. A -3 slope 
is expected for the equilibrium region in both deep and shallow water 
and is consistent with the wave slopes. 
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Figure 10. Sea Surface Elevation spectra at Wave Gauge #1 . 
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Figure 11. Sea Surface Elevation spectra at Wave Gauge. #3 
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V. CONCLUSIONS 



The values of the pdf's of sea surface elevation were positively 
skewed verifying the asymmetrical shape of the waves in the breaker zone. 

On the other hand, the skewness of the velocity components were both 
positive and negative and presented no readily apparent pattern. 

The mean values of the onshore and alongshore components of flow indi- 
cated the presence of a nearshore circulation cell probably associated with 
the rip current which was frequently observed to the south of the instru- 
ments. 

The values of the velocity energy-densi ty spectral components calculated 
from wave spectra using linear theory indicate a qualitative, but not 
quantitative, relationship. In general, linear theory overestimated the 
magnitude of the velocity spectra because of wave directional ity. 

The high coherence between waves and onshore flow indicated that most 
of the motion in the breaking wave is wave-induced. The coherence between 
waves and onshore flow was used to separate the turbulence to wave-induced 
velocity components. The measurements showed that over the range of col- 
lapsing to spilling breakers a reasonable value for the ratio of turbulent 
to wave-induced velocity is approximately 0.75. 

The waves lead the onshore flow in spectral phase on the average 
by less than 20 degrees, implying that the "curling" crest of the wave 
arrives prior to maximum water particle velocity. 

The wave and velocity energy-densi ty spectra at higher frequencies 
exhibited slopes of -5 and -3, respectively. This is indicative of deep 
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water wave conditions. A possible explanation for the surprising spectral 
slope is that most of the waves measured were either spilling breakers or 



reformed waves, 
tion conditions 



It may be that the measured waves did not attain satura- 
during breaking, as do plunging and collapsing breakers. 
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DEPTH OF IMMERSION (FEET) 



APPENDIX C 
CALIBRATION FACTORS 



DATE 


INSTRUMENT 


CALA* 


CALM** 


MARCH 
16,21 ,23 


WG #1 


-1.154 


.00474 


17,19 


WG #1 


-1.284 


.00503 


18 


WG #1 


-1 .453 


.00520 


ALL 


FM #1 






DATES 


ONSHORE 


.009 


.00764 




ALONGSHORE 


- .006 


.00747 




FM #2 








ONSHORE 


+ .131 


-.00911 




ALONGSHORE 


- .084 


-.00818 




WG #3 


- .2798 


.00156 




FM #3 








ONSHORE 


.007 


-.00820 




ALONGSHORE 


.006 


.00763 



*Calibration Additive Factor for wave heights is given in meters; 
CALA for flow is given in meters per second. 

**Cal ibration Multiplication Factor for wave heights is given in 
meters per bit; CALM for flows is given in meters per second per bit. 
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